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Abstract 

Multi-loop superstring amplitude are calculated in the conventional gauge where Grassmann 
moduli are carried by the 2D gravitino field. Generally, instead of the modular symmetry, the 
amplitudes hold the symmetry under modular transformations added by relevant transforma- 
CN | tions of the 2D local supersymmetry. If a number of loops is larger than 3, the integration 

measures are not modular forms. In this case the expression for the amplitude contains an 
, Q ■ integral over the bound of the fundamental region of the modular group. 

1 Introduction 

In the Ramond-Neveu-Schwarz theory the world sheet is often specified [1] as the Riemann 
surface with spin structure |2j. The spin structures are not invariant under transformations of 
the 2D supersymmetry. It leads to the well known difficulties [IJ [3j H] in the calculation of the 
multi-loop interaction amplitudes. They depend [1] on the 2D gravitino field [3j 0]. It means 
that the world-sheet supersymmetry is lost. Indeed, in the superstring theory the vierbein and 
1 the 2D-gravitino field are the gauge fields of the group of local symmetries on the string world 

sheet. Owing to gauge invariance, the "true" amplitudes are independent of the choice of the 
(Sj ■ gauge of the above fields. 

In j5] true two-loop amplitudes have been obtained. The calculation of the Ramond-Neveu- 
Schwarz amplitudes with any number of loops has been done [61 [7] in the supercovariant scheme 
[8j [9l[T0]. In this case the zweibein and 2D gravitino field are conformally flat. The string world 
sheet is specified as the (1|1) complex, non-split supermanifold. The supermanifold carries a 
"superspin" structure 0, QUI E] instead of the spin structure [2] . The superstring amplitude is 
obtained by a summation over the superspin structures. The superspin structures are super- 
symmetric extensions of the spin structures [2]. In this case the twist about (A, £?)-cycles is, 
generally, accompanied by a supersymmetric transformation including fermion-boson mixing. 
The fermion-boson mixing arises due to the presence of Grassmann moduli that are assigned to 
the (1|1) complex, non-split supermanifold in addition to the Riemann ones. The fermion-boson 
mixing differentiates the superspin structures from the ordinary spin ones. Indeed, the ordinary 
spin structures [2] imply that boson fields are single-valued on Riemann surfaces. Only fermion 
fields being twisted about (A, 5)-cycles, may receive the sign. The g-loop spinning (fermion) 
string interaction amplitude (with g > 1) is given by an integral over (3g — 3\2g — 2) complex 
moduli and over interaction vertex coordinates on the supermanifold. The integrand (the local 
amplitude) has been explicitly calculated [6] for every superspin structure. The calculation em- 
ployes only the gauge symmetry of the fermion string. In doing so, the partition functions are 
computed from equations [H [12] that are nothing else than Ward identities. These equations 
realize the requirement that the superstring amplitudes are independent of both the vierbein 
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and the gravitino field. Therefore, the obtained multiloop amplitudes are consistent with the 
gauge invariance of the superstring theory. The world-sheet gauge group is so large that the 
Ward identities fix the partition function up to a constant factor. This factor is determined by 
the factorization condition following from the unitarity equations. The module space integral 
of the local amplitude is, however, ambiguous [HE] under non-split replacements of the moduli. 
The ambiguity is also present in the sum over super-structures for the superstring amplitude. 
The ambiguity in the superstring amplitude is resolved j7] so that the cosmological constant is 
equal to zero. The obtained superstring amplitudes are finite. The one-, two- and three-point 
massless boson amplitudes vanish in accordance with the space-time supersymmetry. 

In this paper, we consider the calculation of the same multi-loop amplitudes in the above- 
mentioned gauge pQ where the (1|1) complex supermanifold is split, in the sense that fermions 
are not mixed to bosons under twists around non-contractible cycles. The genus-g supermani- 
fold is specified as the genus-g Riemann surface £ 9 with the given spin structure. The spinning 
string amplitudes each is represented by an integral of a local amplitude where the integration is 
performed over vertex coordinates and over the moduli. The superstring amplitude is obtained 
by summing over spin structures. Grassmann moduli are carried by the 2D-gravitino field <p m 
which is usually specified such that 7 m m = where 7 m is the Dirac 2D-matrix. Non-zero 
components <fi±(z, z) of the (fi m field are given by 

2g-2 2g-2 
4>-{?i Z) = E ^s(j)s-(z, Z) , (j)+(z, Z) = ^s<Ps+{.Z, Z) , (1) 

8=1 S=l 

A s and A s being Grassmann moduli, and (p ST (z, z) fields may depend on the Riemann moduli 
too. As far as the world-sheet supermanifold is split, one seemingly avoids the ambiguities [HIT] 
complicating the calculation of the amplitudes in the superconformal gauge. Also, it seems 
that the amplitudes might possess the modular symmetry and be represented through theta- 
functions and modular forms. Indeed, in the two-loop calculation (5] the genus-2 integration 
measures are modular forms, and the GSO projections of the local amplitudes with less than 
four legs are equal to zero. The GSO projection of the four-point local amplitude does not 
depend on (fi ST . The spinning string amplitude ceases to depend on (f> ST due to the integration 
over vertex coordinates. The papers [5] have initiated the efforts [13J to build genus g > 2 
amplitudes assuming certain properties of the amplitudes, the modular symmetry being among 
them. This strategy meets with difficulties [141 [15], at least for g > 3. 

The calculation of the multi-loop interaction amplitudes in the present paper is similar to 
the calculation [6] in the supersymmetric gauge. It exploits the gauge symmetry on the string 
world sheet and uses no assumptions. In this case the spinning string amplitudes are inde- 
pendent of local variations of the <p ST fields, but the local amplitudes depend on (j) ST (the last 
takes place even in the two- loop case [5]). Integration of the local amplitude over the moduli 
and over vertex coordinates is performed at fixed <ft ST . The integral must be invariant under 
re-definitions of the non-contractable cycles on the string world sheet. The re- definitions of 
the non-contractable cycles are accomplished by modular transformations, but these transfor- 
mations, generally, change <p ST . Returning back to the original S=F is achieved by an extra 
transformation of a local 2D supersymmetry. Therefore, the symmetry group of the amplitude 
consists of modular transformations accompanied by the relevant supersymmetry ones. These 
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supermodular transformations are conveniently discussed in the supersymmetric description 
[12J of the fermion string on the (1|1) complex supermanifold. The period matrix [6j [7J [11] 
on the above-mentioned supermanifold collects periods of scalar superfunctions which vanish 
under the supercovariant Laplacian [7J . The supermodular transformation changes |H El EEE] 
this matrix just as the relevant modular transformation changes the period matrix [TH] on the 
Riemann surface E g . 

A loss of the supersymmetry in [TJ occurs because the difference between the supermodular 
and modular symmetries was ignored and, also, because of an incomplete calculation of the 
ghost zero mode contribution to the integration measure, see Sec. 2 and Sec. 3 below. 

Since the superscalar functions depend on m , the period matrix on the (1|1) complex 
supermanifold is, generally, distinguished from the period matrix on E 5 by terms proportional to 
the Grassmann moduli. In this case the integration of the local amplitude over the fundamental 
region of the modular group leads to the loss of 2D supersymmetry. As the result, the spinning 
string amplitudes depend on <p ST . To restore the supersymmetry, the integration over the 
fundamental region of the modular group must be supplemented [7J by the integral around the 
boundary of the region. If g < 3, the periods of the superscalar functions can be taken [5] as 
the moduli set. In this case the boundary integral does not arise. If g < 3 and the moduli 
setting (TJ is used, the boundary integral is removed (see Sec. 4) by a re-definition of the local 
amplitude. The integration measures are given by modular forms for both g = 2 and g = 3. 
Unlike the two-loop case, the GSO projection of the four-point, three-loop amplitudes ceases 
to depend on T due to the integration over vertex coordinates, just as it arises in each of 
the spin structures. If g > 3, periods of superscalar functions depend on Grassmann moduli 
for any choice of moduli variables. The boundary integral is present in the expression for the 
amplitude, and the integration measures are not modular forms. It is akin to what occurs in 
the superconformal gauge. Hence the strategy [TJ] is not in accord with the 2D supersymmetry. 

In Sec. 2 the integration over moduli is discussed. In Sec. 3 local amplitudes are calculated. 
In Sec. 4 two- and three-loop amplitudes are considered in more details. 



2 Integration of local amplitudes 



As noted in the Introduction, the period matrix on the (1|1) complex supermanifold determines 
the periods of the scalar superfunctions. The scalar superfunctions vanish under the super- 

depend on the gravitino field (IT 



Laplacian (D^' 1 — D^ T D+)/2 where T denotes transposing. Operators and D_ 



We assume that the S=F fields do not overlap. Then [7j 



D 



(</>) 



D + -</>+(z,z) 



D 



(<t>) 



4 + W 3E 



D = dd z + a. 



(2) 



where •§ is the superpartner of z. The scalar superfunctions J'( B/ *\z,z,'&) are associated with 
the right movers, and (z, z, *&) are associated with the left movers; a and a' mark spin 
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structures. The desired sup erf unctions are represented as 

J^\z,z,$) = 4 R °\z,z)+$4 R(T \z,z) , J^'\z,z,€) = 4 L °'\z,z)+$ v l LcT '\z,z) . (3) 
They can be found from the equations 

D ( f J™ (z, z,#) = 0, D^jM (z, z,&) = 0. (4) 



lf(f) + (z, z) = <ft-{z,z) = 0, then (z) is reduced to the scalar function J r (z, z), and .1^^ ( 



z,z) 



is reduced tqj J r (z). Under 27r-twists about 5-cycles on the Riemann surface E 9 (which specifies 
the genus-g supermanifold in question) the J^'^z^z) functions receive periods forming the 
Qyto) matrix. Correspondingly, the periods of j( La ^(z,z) form the Q( LtT ) matrix. In this case 

where Q( Ra ) and Vt <yLcr ' S) vanish when all the Grassmann moduli are equal to zero, and Vt is the 
period matrix on H g . Eqs. @ can be transformed to the integral equations. In doing so, the 
desired equations for the (j( Ra \z, z), 4 {z, z)) pair and for Vt^^ elements of -matrix 
are found to be 

J^(z, z) = J r (z) - ±- [ d z , \n[E(z, z')) <P_{z', z')4^\z', z')d 2 z' + const , 
4^\z,z) = -L I S a (z,z')<P4z',z>)d z ,4 R °\z',zi)d 2 z, 

Z7T J 

= Q sr - i [ d z J r (z) <f>_(z, z)4^\z, z)d 2 z (6) 



where J r (z) is the scalar function on E 9 having the Q nr periods, E(z, z') is the prime form and 
S,j(z, z') is the Szego kernel [U [T7]. For the even spin structure a = (<7i, 02) it is given by|TJ 

, 9[a](z — z') , f z ' 

S a {z,z') = — , z-z = v(x)dx, v = {v s (x)}, v s (x) = d x J s (x) , 



E(z }Z i)6[a}(0)> 

9[a](z) = 6{z + fia"! + er 2 ) exp[i7ro"if2(Ji + 27rio"i(z + cr 2 )] , o- x = {cr is }, 1 < s < g (7) 

where z is related to z C S 9 by the Jacobi mapping. Further on, #[cr](z) is [U [TF] the theta- 
function with characteristics = (01,02) corresponding to the spin structure 0. The 6[a]{z) 
function is related with the Riemann theta function 6(z) as it is seen from the second line 
of (J7J. The first two equations in are equivalent to the first equation in (J4j). To verify 
it, every integral equation reduced to the differential equation by the d z operator. Indeed, 
S a (z,zj -»■ 1/(2 - z% d z ,\n[E(z,z')} -> -l/(z - z') at z ->■ z', and = ir5 2 (z) where 

5 2 (z) = 5(Re z)5(Im z). Therefore, the integration equations (jl]) are reduced to the differential 
ones. Due to (j2]), these differential equations are identical to the first equation in (J3J. The third 
equation in ([6]) determines the periods of J^'iz^z). It is obtained from the first equation by 
means of the transformation which is assigned to the 27r-twist about the B s -cyc\e on E 5 . Since 
the kernels in (jH]) are proportional to the Grassmann moduli, eqs. (jH]) are solved by the iteration 



1 Throughout the paper, the line over denotes complex conjugate 



4 



procedure. The functions and the period matrix for the left movers are calculated in the similar 
manner. 

Thus £7^°") and QS L(T ' depend on Grassmann moduli. Under these conditions, the moduli 
space integral over the fundamental region of the modular group [TB] is not invariant under the 
supermodular transformations that leads loss of the 2D supersymmetry. To restore the super- 
modular symmetry, the discussed integral is supplemented by an integral over the boundary of 
the integration region. To derive this boundary integral, it is useful to define a function which 
is an extension of the step function p(x) (being p(x) = 1 for x > 0, and p(x) = for x < 0) 
to the case when x = x& + x s contains the "soul" part x s that is the part proportional to the 
Grassmann parameters. Then p(x) is understood in the sense that it is the Taylor series in 
x s . In the calculation of the Taylor series one employs the known relation dp(xb)/dxb = S(xb) 
where 5(x) is the Dirac delta-function, and the "body" Xj, of x contains no Grassmann param- 
eters. Under this convention the fermion string interaction amplitude A a y can be represented 
as the integral of the A a y local amplitude as follows 



A a y = J A^O{n^M L ^) (j[0(z t , Zi ;q,q)d 2 z^ d 2 qd' 



A 



where A a ^< depends on the Riemann moduli (q,q), the Grassmann moduli (A, A) and on the 
(zi, Zi) coordinates of the z-th interaction vertex. In this case q = {q n } and A = {A^}. Generally, 
= fi( Rff )(g,A) and fi (La ' } = ^ La '\q, A). The 6j{z h Zj\q,q) factor is a step function 
product restricting the integration region (that is the fundamental region of the Klein group) 
on the complex ^-planes. The integration region over the moduli space is "restricted" by the 
0{Vt^ R<T \ Vt^ Lr7 ') step function product. In this case 

o(n^\n^) = HpiGj) , Qj = g^isi^M^) . (9) 
j 

The set of the Qi(Q, Q) = conditions gives the boundary [16] of the fundamental region of the 
modular group. The step functions p(Gj) in are treated as the Taylor series in Q p R ^ and in 
matrix elements. Therefore, 

oisi^M L ^) = B Ra B La ,o(n,U) (io) 

where the differential operators 3^ and B^i are defined as follows 

Eu = i + \^& Ra) -^— + - \^\^n {Ra) n {Ra) — — + 

DRa l ^ 2^ lL pq an + nZ^Z^ 1 V9 "rs oq oq "T • • • > 
p<q u ^ L pq ^ p<qr<s U\L pq u±L rs 

p<q VlL pq p<qr<s UZL pq UlL rs 

The derivatives in (|TT|) are calculated assuming that the period matrix elements are unrelated 
to each other up to the transposing operation. Under the operators (flT!) . the step function 
in (flOl) receives <5-function-type terms that leads to the appearances of the integral over the 
boundary of the fundamental region of the modular group. 
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Under the change of integration variables in (jSj), the arguments of the step functions are 
correspondingly replaced. As the result the amplitude (JS} is independent of the choice of the 
integration variables. Strictly speaking, the last statement implies that the integral dSJ) is 
properly regularized at the points where the Riemann surface is degenerate, but we do not 
discuss this matter in the present paper. The amplitude (jSJ) can be derived [7] by a change 
of integration variables in the expression for the same amplitude in the superconformal gauge 
[U [9]. Hence (jHJ) is independent of (f> ST (we have directly verified it for g < 3). 

3 Local amplitude A a , a f 

To derive A a ,a' in ©, we start [SI 13 CE2] with the integral [18] over all the fields, including 
the zweibein and the world-sheet gravitino field. The integral is divided by the volume of the 
local group Q of the world-sheet symmetries of the fermion string. So far as the zweibein and 
the world-sheet gravitino field are arbitrary, we can map [T2l E] the Riemann surface onto the 
complex plane w choosing the same transition group Qt for all surfaces of the given genus-n. 
There is no integration over any moduli. The zweibein and gravitino fields can be reduced to 
the full set of the reference fields. It is performed by globally defined transformations of the 
Q group that do not change the Qt transition group. The reduction is impossible within the 
full set of the reference fields. The zweibein and the gravitino field are represented in terms of 
the reference fields and of the gauge functions. Since the gauge functions correspond to the Q 
group transformations, the Q t transition group is unchanged and, therefore, it is the same for 
all the genus-g surfaces. In this case the reference fields (for g > 1) depend on {3g — 3\2g — 2) 
complex moduli (defined up to the supermodular transformations). Locally, the reference fields 
are arbitrary. The integration over the zweibein and the gravitino field is transformed to 
the integration over the gauge functions and the moduli. In doing so the Jacobian of the 
transformation is represented by the integral over the ghost fields and over (3g — 3\2g — 2) 
global complex variables dual to the (3g — 3\2g — 2) complex moduli. Calculating alterations 
of the integral under infinitesimal local variations of the reference fields, one can derive the 
Ward identities [SI E] from the condition that the amplitude flS]) is unchanged under the above- 
mentioned variations of the reference fields. The obtained Ward identities are transformed to 
the desirable gauge of the reference fields. The Ward identities can be used for the calculation of 
the local amplitude. Indeed, the direct calculation the amplitude from the integral over the fields 
is hampered as determinants of the differential operators appear in the calculation. Therefore, 
the integral requires a regular izat ion ensuring the independence of the amplitude (jSJ) from 
infinitesimal local variations of the reference fields. In the considered gauge pQ the discussed 
uncertainty is, however, appears in the local amplitude as the factor which is independent of the 
Grassmann moduli and of the vertex coordinates. Excepting this factor, the local amplitude 
can be obtained from the discussed integral over the fields and over the (3g — 3\2g — 2) global 
complex variables (A^JA{) dual to the (3g — 3\2g — 2) complex moduli. The integral is, as 

follows [SI m E2] 

= A { :1> < I] Vi > <t> = J (DJ^)d 2 AW (V) exp[S m + S* + S% + S^ h + S%\ (12) 
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where A^\, is a "vacuum local amplitude" and < Oj Vj ><t> is the vacuum expectation of the 
vertex product (V). The vacuum expectation is calculated in the gravitino field ([TJ). Further, 
(DJ 7 ) is the product of differentials of the fields, the fields being the ghost complex fields and 10 
scalar x N fields with their superpartners (ip N , 4> )• The ghost complex fields are (2, — l)-tensor 
fields (6, c) and (3/2, — l/2)-tensor fields (j3, 7). Further on, 

2 



S m = - V] [ d 2 z 

71 N J 



—dx dx N + r ip dilj N +'ilj dip N + (j)^ip dx N + (p + ip dx N 



qR 
dgh 



- I d 2 z 



71 



-bdc + (3d-f - -0_(&7 + (3dc) + /3(<90_)c 



- f d 2 z -J2[bd<; m + -cj>-^ m - m-Sm- P-^-)A b m + J2^s-M 



(13) 

(14) 
(15) 



where df = d z f and df = c^/ for any function / = f(z,z). Also, </> s _ = (f) s -(z,z) and 
0_ = <j)_(z,z), see ([1]). The and in ffl2|) are obtained by the complex conjugation of 
(fT4|) and of (ITS"]) together with the 4>-(z,z) —> (p + (z,z) replacement. 

Eq. (fl5l) contains a function q m = q m (z,z) that is one- valued under rounds about A-cycles 
and has a discontinuity under twists about S-cycles on the Riemann surface. Let assign the 
z — > g(z) replacemento to the 2tt twist about 5 s -cycle. Then 



^m{9s{z),g a (z)) = S m (z,z) 



dz 



+ 



dg s {z) 
dq m 



(16) 



The last term on the right side of ( !T6|) is the discontinuity. Due to the discontinuity in s m (z, z), 
the integration over zero modes of b- fields in (Tl2l) is convergent. As explained below, the result 
of the integration in (IT2"|) does not depend on a further specification of q m (z,z). 

The V vertex in (fT2l is built using the supercovariant operators ($2§, but it can be verified 
that the T -dependent terms in (T5]) do not contribute to the amplitude. It appears due to 
motion equations following from ffl3|) . So the conventional vertex (19] can be employed. 

With the exception of the (A s , X s ) independent factor in A^ a ,, the amplitude ( TT2|) is ex- 
pressed in terms of the correlation functions. They are calculated from the integral ([12"]) at 
vanishing (A S ,A S ). In so doing the linear sources of the fields and of the global variables are 
added to the exponent (it is the known trick in the calculation of correlation functions). From 
ffl~4]) and f fl5|) . it follows that c = c(z, z) and s m (z, z) are combined into the c(z, z) field where 
c(z, z) = c(z, z) + J2 S S m ( z , z)^m- The < cb > correlator (at X s = 0) is (6] [7] ITT] 



where G b (z,z') — > l/(z 

G b (g s (z),z') 



< c(z, z)b(z', z') >= G b (z, z') 
z') at z — > z', and 



(17) 



dz 



(3(9-1) 

-G b (z,z')+ £ 

m=l 
Xm(z) = ~ 



. , ~Xm(z ) , 

< b(z, z)A b v > 



(19) 



2 For the sake of simplicity we assume the Schottky description of the Klein group, but the using of the 
Schottky moduli is not implied. 
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The last term on the right side of the first equation in ({18]) appears due to the discontinuity (JT6]) 
of <, m . Furthermore, Gb(z, z') is not changed under rounds about A s -cycles. As the function of 
z', the correlator (fT7|) is the conform 2-rank tensor. These properties are sufficient to determine 
both Gb(z,z') and 2-rank-tensor zero modes Xm( z ). Henceforth the amplitude ( 1T2"]) does not 
depend on details of q m (z,z). It should be noted that Gb(z,z') is given [7J [TT] by a Poincare 
series that is not expressed through a local combination of theta-like functions. 

In p] the correlator f lTT]) is mistakenly replaced by < c(z, z)b(z', z') >= Gb(z, z'\p) depend- 
ing on 3(g — 1) arbitrary points p = {p m }. The Gb(z,z';p) correlator is one-valued on the 
Riemann surface, has poles at z = p m and vanishes at z' = p m . This properties determine 
Gb{z,z';p) up no a numerical factor. In our normalization of the fields the correlators are 
related, as follows 

G b (z, z'] p) = G b (z, z') - G b (z, p a )Xm(z') , Xm{Pn) = S mn (20) 

m 

where Xm( z ) are the 2-rank tensor zero modes which are normalized as it shown in (120]) . Due 
to eq. (1181) . the right side of (120]) is one- valued on the Riemann surface. In addition, it has poles 
at z — p a , vanishes at z' = p a and goes to l/(z — z') at z — > z'. Thus the right side of fl20|) 
coincides with G (z, z';p). 

To clarify discrepancy between the amplitude ( Tl2l) and the corresponding amplitude in P 
we transform (fl2l) to an integral where the (6, b) fields each vanishes in 3(g — 1) points on the 
Riemann surface. For this aim the (? m , q m ) functions are properly specified, and the proportional 
to Ph. b m and to (3A m terms are removed from the exponent in ffl2l by a relevant shift of the 
(7,7) fields. In more details, 7-^7 + Em^»A^ where q m = q m (z,z;p) and q m = q m (z,z;p) 
depend on the p = {p a } set of 3(g — 1) points p a on the Riemann surface. Furthermore, 

q m (z,z;p) = / G a (z, z'; {(j)})[d z/ (f)(z' , z') q m (z, z;p) - -cj)(z' ', z')d z ,<; m (z, z; p) 

+d qm (t)-{z,z)}d 2 z', 

q m (z,z;p) = Y\G h (z,p a )N am + — / G b (z,z')cj)(z', z')q m (z,z;p)d 2 z' , 

■Nam ^ t N an Nnm i N na Xn iPa ) ) 
n 

Mum + ^~ [ Xn{z')<P{z', z')Zm(z, Z] p)d 2 z' = 1 (21) 

where A^^ 1 is the element of the iV -1 matrix inversed to the matrix N. The 2-rank-tensor zero 
modes Xn(z) are the same as in ([19]) . The G a (z, z'\ {0}) function satisfies to the equation as 
follows 

d- z G a {z, z; {4>}) = n5 2 (z - z') - ^ <t>- m (z, z)x f r (z') , / x f r {z)<j)_ n {z, z)d 2 z = 5 mn (22) 

m J 

where 5 mn is the Kronecker symbol and x((z') is 3/2-rank-tensor zero modes. In this case 
G a (z, z'\ {</>}) = — < -y(z, z) f3(z', z') > where < -y/3 > is the correlator at vanishing Grassmann 
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moduli. Due to the ~ A{ terms in (fT5|) . it is a functional of S _. The (7, /3)-dependent part 
of the exponent is represented in the kindred manner. The (7, (3) and (7, (3) correlators at 
arbitrary Tm are expressed through the correlators in the case when 0_ m and cj) s+ each is 
localized at z = r s and, respectively, at z = r' s , as follows 



'z, z) = 5 (z - r. 



'z,z) = 6 2 (z-r': 



(23) 



In this case we denote the (7,/?) correlator as < j(z, z) P(z', z') >= —G a (z, z';r) where r = 
{r s }. This correlator was calculated in pp. It has the pole at z = 77, the residue being the 
3/2-rank-tensor zero mode x{( z ') satisfying to the condition x{( r j) — — $sj- Then 

G a (z,z';{0}) = G a (z,z';r) + J2 - f G (T (z,z 1 ;r) ( p s 4z 1 ,z 1 )d 2 z 1 x f Az') , (24) 

x{{z) + -Y.lx{(z')<P r .(z\z')d 2 z'xi(z) = Q. (25) 

Indeed, one can verify that the right side of (125]) satisfies to (1221) . The kindred expression exists 
for the (7, (3) correlators. Once the integration over the global variables being performed, the 
vacuum amplitude A^\, is found to be as follows 



A 



(v) 



{DJF) W R(T (p, p) W La , {p', p') exp[S n 



where 



detTV" 



det [Xm(Pa)\ 



3(9-1) 

II b (Pa,Pa) 
a=l 



r 2(g-i) 

n * 



7T 



P(z,z)</>j 



iz, z 



d 2 z 



(26) 



(27) 



and Wl<t>{p' ,p') is the kindred expression associated with the left movers. Elements M n m of 
the M matrix is defined in (T2TT) . Zero modes Xm{Pa) are defined by f lT§|) . The product over 
a in Wn a (p, p) provides the vanishing of the b(z, z) field at z = p a . Eq. f )2T|) differs from the 
corresponding expression in [TJ by the det M factor (apart from the fact that in [1] the set of 
2-rank tensor modes is not specified). This factor generates the terms in Wr„(p, p) which are 
added to the vacuum expectations of the supercurrent products (arising from the expanding 
of exp[S^] in As the result, the Gb(z,z';p) correlator is replaced by the Gb(z,z') one. If 
0_ depends on the Riemann moduli, additional terms in the amplitude also appear due to the 
<9 gm 0_ terms in f)15p . The kindred thing arises in the integral over the (7, (3) fields. 

As was noted, the vacuum local amplitude f|T2|) contains an uncertain factor. The factor is 
independent of the (A s , X s ) moduli. Nevertheless, it depends on the 4> ST fields ([1]) because 
and contain the derivatives of the gravitino field with respect to the Grassmann moduli. The 
factor can be calculated from Ward identities [7]) as was mentioned above. For the commonly 
used setting fT25|) for (f> ST the discussed factor was already calculated in [I]. If this setting is 
employed, the local vacuum amplitude is represented as follows 



A 



(«) 

a,a f 



1 5 



det ImVt 



Zriq, r)Z a ,(q, r')Z ( ^f (q, q, A, A, r, f')Z^\q, A, r)Z { f\q, A, r>) {21 
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where r = {r s }, r' = {r' s }. The (A, A)-independent factor Z a (q,r) can be taken from pQ. The 
other three factors differ from the unity only because of proportional to (A, A) terms. Among 
of them, Z^ T f t \q, q, A, A, r, f') is due to the expanding in (A, A) of exp S m . The last factors are 
due to the (A, A) expanding of the rest exponential in ( fl2l) . The calculation of Z^ h \q, A,r) is 
ambiguous because, as it was noted above, the correlator < j(z, z) 0(z', z') >= —G a (z, z';r) 
has the pole at z = Vj. To resolve the ambiguity^! one treats fields ( 1231) as the limit of a spread 
fields (p ST using eq. (l2^]) to calculate of the (7,/?) correlator. 



4 Two- and three-loop amplitudes 

In the discussed case the periods Q nm = fl mn of the scalar functions can be taken as moduli. 
By using eq. ({TO]) , the integration by parts is performed in ([8]), and A aa i is represented as 

A^, = J 0(Q, n)B T Ra B T La , L CT v II &(*i> *)] d " X IT d2 *i ( 29 ) 

i i 

where 0(zi, zA = 0(zi, Zi, tt, ft), and (B^, B^,) are obtained by transposing of (Bru, Blct')- If 
the Schottky description is employed, then 

d J (z J ,z j )=p(l-\g'f)p(l-\~g' s \ 2 ) (30) 

where g' s (z) = dg s (z)/dz, the z — > g s (z) transformation is assigned to 27r-twist about -B s -cycle. 
The z —> g s (z) transformation is inverse to the z —> g s (z) one, so that g s (g s { z )) — 9s(g s ( z )) = z. 
As above, p{x) is the step function. Instead of the boundary integral in the moduli space, eq. 
(129]) contains the integral along the boundary of the fundamental region of the Klein group. 
This integral arises due to the action of the (B^, B\ a t) operators on the 0(zi, Zj) step function 
products. This boundary integral can be reduced to the integral over the fundamental region 
of the Klein group. The above-mentioned reduction is performed employing the set of certain 
functions U pq (z) = U qp (z). The U pq (z) function is unchanged under twists about A s -cycles, but 
it has the discontinuity under 2-7r-twist z — > g s {z) about 5 s -cycle, as follows 

U Mz)) = ^U pq (z) + ^. (31) 

The last term on the right side of ( 13T1) is the discontinuity of U pq (z). Then 

A^ 90 ^^ d\ = J A^MM 90 ^ <? Zl . (32) 

Indeed, from (|30|) . the left side and the right side of eq. (|32|) each is the integral along the 
(A s + A s ) contour where the A s contour is given by the ^^(z)! 2 = 1 condition, and the A s 
contour is given by the |(^(z)| 2 = 1 condition. By the z — > g(z) replacement the integration 



3 In pQ this matter is treated inexactly. 
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along A s is reduced to the integration along A s . The right side of (152]) is calculated using 
eq. (fT8j) and taking into account that A a y is (1, l)-tensor in (z^ Zi). The left side of ( 132|) is 
calculated using the relation 

<9P(1-I^)| 2 ) _ <9p(l-|^)| 2 ) , dp{l - \g' 8 (z)\*) ^ 



Lpq 



z^ 9 s(z) dVt pq 9' s (z)dz 



As the result, the same expression appears for both the left and right sides of (1521) that proves 
the validity of eq.(J32J). The integration by parts being performed, eq. (l2"9|) is reduced to the 
integral over the fundamental regions of the modular group and of the Klein group (there is no 
any boundary integral), as follows 

4t,<t' = / n d 2 n mn n d 2 Zj , = j A { ™, d) n , a^ = bia^, m 

m<n j i 

where A a)(T ' is given by ( IT2|) . Operators B^ and B\ a , are obtained by the d/dVt pq — > V pq 
replacement in £>^ CT and, respectively, by the d/dVL pq — >■ P pg replacement in B\ a ,. In this case 

p « = A; q ~ ? • (35) 

The U pq (z) function is constructed using Gb(z,z') and zero modes Xnl- 2 ); see P^jh ([HD and 
( TT9]) . In g = 2 or g = 3 cases n- index is replaced by a pair (jl) of indices listing the period 
matrix elements. Then we use notation X(ji)( z ) instead of Xn( z )- As it is shown below, 

X(jr)(z) = -2mVj{z)v r {z) (36) 

where v s (z) = d z J s (z) is 1-form. For U pq {z) one can use a sum of Gb(z,Wj)N^ pq ^ over 3(g — 1) 
arbitrary points Wj where iV -1 is the inverse to the N matrix (detiV ^ 0), whose matrix 
elements are N( pq )j = X(pq)( w j)- This U. pq {z) depends on 3(g — 1) arbitrary points Wj. 

For the calculation of GSO-projection it is more convenient to use alternative functions which 
depend on {2g — 3) points Wj (these points can be identified with the 2D-gravitino location 
points). These functions are built using 2-rank-tensor modes (j(u), ( s (u) and t{u) which are 

Q( z ) = -Y,r' j {w s )Uz) - [r>0 -gr;K)^K)]f(z), 

s=l s=l 

Uz) = t s (z)-t' s '( Wi )t(z) (37) 

where 2-rank-tensor modes Tj(z) and f s (z) satisfy conditions that Tj(wi) = 8ji, f s (wi) for all 
wi while t' s (wi) = 5 s i for I = 1 and (if g — 3) for 1 = 2. In this case, for any f(x) function, 
f'(x) = d x f(x) and f"(x) = d 2 f(x). For g = 2 one can set Ti(z), f\{z) and f(z), as follows 

vf{wi) vi{wijd'{wi,wij 2 \d'(wi,wi) / 
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where d'(x,y) = d x d(x,y) and d(x,y) = Vi(x)v2(y) — Vi(y)v2(x). For g = 3 one can set 

d(z,w 2 )d(z,w 3 ) d(z,w 2 )d(z, w 1 ,w 2 ) 

T i\ z ) = -r, T77 7 > n{z) ~ 



d(wx,w 2 )d(w 1 ,w 3 ) ' " d(w 1 ,w 2 )d'(w 1 ,w 1 ,w 2 ) ' 

where d(x 1 ,x 2 ,x 3 ) = detVj(x{), d'(x,y,z) = d x d(x,y,z) and d(x,y) = Vi(x)v 2 {y) -Vi(y)v 2 (x). 
Functions t 2 (z), t 3 (z) and f 2 (z) are obtained by a replacement of indices in T\{z) and in f\{z). 
In this case the U. nm {z) functions are defined by 

2g-3 g-1 

{z,w l )r nm (40) 

j=l s=l 

where (j,nm, (s,nm and f nm are coefficients of expansion of (j(u), ( s (u) and of r{u) in the 
X(nr)(z) / (— 2m) modes. Derivatives with respect to the Riemann moduli of the correlators and 
the functions are calculated in line with [SI Ej. As is proved below, 

^ dJ s (z) dR(z,w) 

}^ X{mn){w)-^ = -x V a {W) - V s (z)& b {z,W) + /J, S {W) (41) 

where R(z,w) is the scalar, holomorphic Green function [S] for the d z d z /n operator. This 
R(z,w) is not changed under twists about A s -cycles. Under the 27r-twist z — > g s (z) about B s 
cycle it is changed as follows 

R{g s {z), w) = R(z, w) + 2mJ s (w) , R(z, g s {w)) = R(z, w) + 2mJ s (z) . (42) 

The explicit form of /a s (w) is not employed in this paper. The Green function R(z,w) differs 
from the usual Green function In E(z, w) in (j6]) only by the scalar zero mode contribution. Thus 

d z d w R(z, w) = d z d w In E(z, w) . (43) 

Differentiating ( 14T]) . one obtains that 

v- , ,dv s {z) d 2 R(z,w) 

L X(mn)(w)— — = v a (w) - d z [v s {z)G b (z,w)\ . (44) 

To derive eq. (l36|) . the z — > g n (z) replacement in (|4"T|) is performed and the relation J s (g n {.z)) = 
J s (z) + Q sn is used, along with eq. ( ITS]) and with the following relations 



dJJz) 



dq v 

fdJ_ 

V dq m ) z dq m ' dq. 
In eq. ( 1451) the dervative is calculated under fixed z 



dUg n (z))\ d9n(z) (45) 



z^ 9a {z) \ dq m ) z dq 

dJ s (g n (z))\ _ dJ s (z) + dn sn 



12 



To prove eq. (jHJ), the difference A s (z, w) of the left and right sides of (Pill) is represented as 

/du 
A s (u,w)d u R(u,z) — (47) 
Zm 

where the integration is performed along the contour surrounding the point u — z. The integral 
(|47p is reduced to the integral along the boundary (13"01 of the integration region. Further, by 
the z — >■ g n {z) replacement, the integral along the \g' n (u)\ 2 = 1 contour is reduced to the 
integral along the contour \g' n {jj n {z)\ 2 = 1. Then, by using eqs. (fl8|) and (146]) . the integral ([47| 
is reduced to the sum of expressions, each being proportional to the integral of d u R(u, z) along 
the \g' n (u)\ 2 = 1 contours. So far as R(u,z) is unchanged under twists about A s -cyc\es, the 
integral vanish that proves eq. (141 p . 

The calculation of the same derivatives of the scalar field correlator < x(zi, zi)x(z 2 , z 2 ) >= 
—X(zi, z~\\ z 2 , z 2 ) is simplified, if it is chosen as follows 



X(zi, z\\ z 2 , z 2 ) = -ReR(z 1 , z 2 ) + ir Y ImJ s {zx) 

^ s,t 



1 



ImQ 



ImJ t {z 2 ). (4* 



St 



This correlator differs from the scalar field correlator in pQ by a scalar zero mode contribution. 
In the calculation of the amplitude both the correlators can be used on equal terms. Further, 
it can be proved that 

^ dR(z,u) dR(z,w) dR(w,u) dR(u,w) dR(w, z) 

Z^X(mn) w Q$i mn dw dw dw dw 

„ , .dR(z,u) _ . .dR(z,u) . . . « — . « , . 

-G b (z, w) — - -G b (u, w) — + fi(w) + (z + u)fi(w) . 49 
dz du 

An explicit form of fi(w) and of fl(w) is not used in the paper. The proof of eq. ( 1491) is similar 
to the proof of eq.flU]). From flUJ) and (j4"9"j) . it follows that 

X{mn)(w) ' ' — = -4:X(z,z;w)X(u,u;w) - X(z,z;u)G b (z,w) 

—X(u,u;z)Gb(z,w), X(z, z;u) = d u X(z, z;u,u) . (50) 

If the relation 

V V, J?\ 

da 



dT 

Yl XmW^t- = U(z) (51) 



takes place for certain T(z), then the derivatives on its left side are expressed through T-L(wj) 
as follows 

Qf 23-3 g-1 

- 27TZ— = Y %( w i)0,rmi + 9 Ws n(w s )Cs,nm + ^^(wi)^ . (52) 

To calculate GSO-projections, the rj points (and the r- ones) in fl23|) are submitted to the 
conditions 

2g-2 

Y r, - 2 A = (53) 

i=l 
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where A denotes the vector of the Riemann constants and is related to by the Jacobi 
mapping . Any (g — 1) points in the {rj} set can be taken at will. It follows from (153]) that 



9-1 



d(r Pl , ...,r Pa ) = det[v s (r P] )} = , v g {r Pi ) = J2 a s v s (r Pi ) 
where a s are the same for every r Pi C {r^}. Furthermore, 



(54) 



s=l 



9-1 



d(r ai ,...,r ag _ 1 ,z) = f(z;g)d gg (r SlJ ...,r Sg _ 1 ), f(z;g) = v g (z) - ^a s v s (z) (55) 



8=1 



where d gg (r Sl , . . . , r s x ) is (gg)-minor of the d(r Sl , . . . , r. _ x , 2) determinant. Under (153"]) . the 
Z a (q,r) factor in (|28|) is simplified essentially. Besides, the < 7(2, z){3(z', z') > correlator can 
be represented as follows 



< j(z, z)P(z>, z') >= -ii^ls a (z, z>) . 

f{z;g) 



(56) 



The sums over spin structures (GSO-projections) are calculated in the known manner using 
Riemann relations and Fay identities (see eq. (45) in Ref. [TZ]). The sums (GSO-projections) 
of the amplitudes vanish, if the vertex number n < 4. If n = 4, then only that part of (V) in 
(|12p contributes to the sum, which contains the product of all the fermion fields. Using PXj) 
and (150]) . one can show that the discussed contribution to A™£® is factorized in the {rj} and 
{r'j} so that -4^./ in (134]) is replaced by an expression 



a (mod) 



1 5 



<y >^)(n,n,A,r,zy,z v )^(fi,n,A,^^,^) (57) 

_det [imiij J 

with zy = {zj}, and < Vq > to be the vacuum average of Vo where 

V = exp[^ z • x(zj, z~j)} , In < V >= - ^ ki ■ kj < x(z iy Zi)x(zj, Zj) > , (58) 

kj being 10-momentum of the j-th boson. Further, terms due to the differentiating of det ImQ 
in ( 128"]) are canceled with that part of Z^,(q, q, A, A, r, f') in ( fl~2]) which arises due to the last 
term on the right side of ( 148]) . As the result, A$(£l, fi, A, r, 2y, z v ) is represented as 



^ 9) (fi,^,A 



1 



<V > 



Bl(Q,X,r,z v ) 



<V >Al 9 \n,\,r,z v ) 



(59) 



where A { J>\n,\, r, zy) is holomorphic in its arguments. The superstring amplitude A g (that is 
GSO-projection of (1341) ) at # = 2 and g = 3 is given by 



det ImVt 



-i 5 _ 

< U > ^> r > ^v, z"v) A?(H ^, r', zv, zv) II ^» II d 2 Zj , (60) 



^(fi,fi,r,^,^) = E I A^\n^\r,z v ,z v )\{dXj. (61) 

/-r «^ 
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The calculation of ( 16T]) is simplified drastically, if one sums up over a before the fl mn derivatives 
will be taken. Under conditions ( )53|) . there are quite simplified the Z a (q,r) factor in ( 128]) . The 
dependence on boson polarizations in ( )6T|) is extracted in the form of the factor \K\ 2 that is 
the same for g = 2 and g = 3, and that is the same as in [5]. 

For the two-loop amplitude the result [5] is reproduced. In this case A g (Q, Q, r, zy, zy) = 
A 2 (Q, r, zy, zy) is found to be 



A 2 (Q,Q,r, Zy,Zy) 



K 



Z 2 (Q,r, Zy) 



Vlia <V > + < dx(n) > v ■ < dx(r 2 ) >y +A 2 



(62) 



where A 2 is independent of boson 10-momenta, x(r,r) = {x (r, r)} and 
Z 2 {n,r,z v ) = 2 ^ _ oW</ _ oN [I f{zi, 2) , < &r(r) >y= ■ 



/'(n;2)/'(r 2 ;2) 



(63) 



In this case fi(r) is 1-form f s (r) for s = 1 and /'(r; 2) = d r f(r;2). The /(r;2) function is 
defined in (155]) . and < Vo<9 r x(r, f) > is the vacuum average of Vo<9 r a;(r, f). The X> operator in 
([52]) is given by 



13 = TEK( r i) w «( r 2) + M r 2)v q {n)] 



p<g 



on 



" - - YMpiW— 



d Zj 



(64) 



where U pq (z) is defined by ( HOI) . The *4. 2 term in (1621) is actually equal to zero. The proof of this 
statement will be presented in a future publication where the summation over spin structures is 
planned to give in details. The other terms in ( l6"2"j) is calculated using (|4"8~]1 . f l50|) . ( 152"]) and f )58|) . 
In this case wi in (138]) is chosen among r 3 -, say, wi = r 1; and the same point W\ = r\ is chosen 
in (J52"j) . Then *4. 2 (fi, fi, r, zy, zy) = *4 2 (fi, zy) in ( 162]) is found to be (below Z 2 = Z 2 (0, r, Zy) 
in eq.(|63D) 



^4 2 (fi, Zy) 



K 

16tt 2 



Z 2 ^ fc» • fcj 



v l {r x ) 



X(z u Zi\ r 1 )X(zj, Zj] ri) - X(^, ^; n)X(z,-, zf, r 2 ) 



(65) 



see (150]) and (158]) for notations. Further, using eqs. (144"]) . to ( HHl) and ( 152]) (for g = 2), one can 
derive some number of identities due to the fact that the left side of (J5"2"j) is independent of u>i. 
In the explicit form the identities will be given in the future publication. Using these identities, 
one can represent (1631) as 



A 2 (£l, Zy) 



167T 2 



16/(^; 2)/(^;2)t;?(r 2 ) 



1 /'(r 2 ;2)0<y o > 

777T - o\ 51 ( 66 ) 



<V >V 4/(^;2) 



where notations are given in (|58p and in f ]6"3"]) . A detailed deriving of (16"B"]) is planned in the 
future publication. The last sum on the right side of (166]) originates in (160]) terms corresponding 
globally defined derivatives with respect to Zj. Thus it does not contribute to the superstring 
amplitude. The first sum on the right side of (166]) is calculated using ( 163]) and the equation 
f'(ri,2)/vf( y r 1 ) = — f'(r 2 ; 2)/v\{r 2 ) which is obtained by differentiating det[u,-(r s )] = with 



15 



respect to r±. In this calculation one employs that dr 2 /dr\ = — vi(r\) / v\{r 2 ) . The last equation 
is obtained by the differentiation of ( )53|) with respect to r\. Once the last sum on the right 
side being omitted and the 10-momentum conservation being taken into account, A 2 (Q,zy) 
becomes to be: 

K K 
A 2 (Q, zy) = k i ' k iVi{z j )v 1 (z l )v 2 {z m )v 2 {z n ) = q^Y({z, z , k Y) , (67) 

Y({z, z, k}) = hy(l, 2; 3, 4) + y(l, 3; 2, 4) + y(l, 4; 3, 2)] , (68) 
6 

y(i,j; l,m) = {h - kj) ■ (ki - k m )d(zi, Zj)d(zi, z m ) (69) 

where d(z^Zj) = Vi(zi)v 2 (zj) — v 1 (zj)v 2 (z i ). The superstring amplitude A 2 is obtained by 
substitution of A 2 (Q, zy) = A 2 (Q, Q, r, zy, Zy) to (!60|) . It coincides with the amplitude in [5]. 

In the three-loop amplitude, A g (Q, Q, r, zy, zy) = As(ft, tt, r, zy, zy) = ^4.3 in (16T|) is calcu- 
lated the kindred manner. Potentially, ^.3 might contain terms of fourth-order, two-order and 
of zero-order in 10-momentum components k^ of the interaction bosons. It can be proved that 
the fourth-order terms disappear in (|6"0"|) due to the integration over Zj. It seems plausible that 
zero-order terms are absent in A3, but it needs a further study. The calculation of quadratic 
in kj' 1 terms in ^3 has analogy with the calculation for g = 2. Unlike the g = 2 case, these 
terms in *4 3 depend on rj though the whole amplitude f )60|) is independent of rj. Besides, the 
discussed terms in A3 contain the Gb function ( TT7l) that, as it was noted above, can not be 
expressed in term of a local combination of theta-like functions. 

The paper is partially supported by Russian State Grant Scientific School (RSGSS)-65751. 
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